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Abstract 

Several aspects of hadron physics are well described by a simple 5D effective field theory. 
Baryons arise in this scenario as "large" (and therefore calculable) 5D skyrmions. We 
extend and refine the existing analysis of this 5D soliton, which is fairly non-trivial 
due to the need of numerical methods. We perform the complete quantization of those 
collective coordinates which are relevant for computing the static observables like the 
nucleon form factors. We compare the result with simple expectations about large- A^c 
QCD and with the experimental data. An agreement within 30% is found. 



1 Introduction and Conclusions 



Certain 5D effective gauge tlieories, often referred to as "Holograpliic QCD" or "AdS/QCD" 
models [1-4], closely resemble low-energy QCD in the limit of large number of colors Nc. 
The similarity is qualitative as these theories contain, like large-A^^c QCD, infinite towers 
of weakly interacting mesons, but also quantitative. Leading order calculations in such 5D 
models typically describe the physics of the lightest mesons to 10% accuracy in terms of an 
extremely limited number of parameters. These results are compatible with the hypothesis 
that leading order calculations in the 5D model reproduce large- iVc QCD. 

Baryons arise in this scenario as solitons with a conserved topological charge which rep- 
resents the baryon number. They are the 5D analog of skyrmions [5,6] (see [7] for a review), 
so we will refer to them as 5D skyrmions. They differ from ordinary skyrmions, however, 
in the important aspect of calculability, as shown in [8,9]. The skyrmion solutions obtained 
in 4D models of mesons -which may or may not contain some vector resonances on top of 
the pion field- have a size which is of the order of the inverse cut-off of the theory, and 
incalculable UV effects do not decouple. This problem, which constitutes the main reason 
of theoretical dissatisfaction about the 4D Skyrme model, is solved in the 5D case because 
the size of the 5D skyrmion is larger than the 5D cut-off. 

A string construction named "Sakai-Sugimoto model" [10] might provide a UV comple- 
tion of the AdS/QCD models and give them an interpretation in the standard framework 
of AdS/CFT. In the low-energy supergravity limit of large 't Hooft coupling A ^ oo the 
Sakai-Sugimoto model reduces indeed to a 5D theory with U{Nf) gauge symmetry {Nf de- 
notes the number of flavors) and two AdSs-like boundaries on which the sources for Left- 
and Right-handed currents are located. This can be rewritten as a U{Nf)L x U{Nf)fi theory 
living on one-half of the space with one AdSs-like (UV) boundary on which both sources 
live and one IR boundary on which symmetry-breaking conditions as in eq. ^ are imposed. 
The Sakai-Sugimoto model is almost equivalent, in the limit in which practical calculations 
are performed, to the effective theory considered in the present paper. There is however 
a difference which, as remarked in [8,9], becomes extremely relevant in the baryon sector. 
In the Sakai-Sugimoto model the effective 5D interaction scale M5 is proportional to the 't 
Hooft coupling A, i.e. M5 oc A, while the coefficient of the Chern-Simons (CS) term, which 
is fixed by the Adler-Bardeen anomaly, has no A factor. This implies, given the definition 
in eq. ([7|,70cl/A— s>0. The parameter 7 controls the size of the skyrmion, p oc 7^/^ in 
the Sakai-sugimoto model and for this reason the string effects, which are encoded in higher 
dimensional operators, do not decouple and there is no advantage with respect to the usual 
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4D Skyrme model in what concerns calculability. It is interesting, nevertheless, to forget 
about higher dimensional operators and study baryons in this framework, describing them 
in terms of "small" Yang-Mills instantons [11] or in terms of 4D Skyrmions [12]. 

In the present paper we complete and refine the analysis of the 5D skyrmions presented 
in [9], with the aim of computing nucleon static observables and in particular the current 
form factors. To this end we need to perform a complete quantization of the 5D skyrmion 
collective coordinates, which is a non-trivial task as it requires to solve numerically a new 
set of partial differential equations. Obtaining predictions for the form factors at non-zero 
transfer momentum requires, moreover, an increased precision of the solution, which we 
obtain by refining our numerical method. The nucleon form factors in the Sakai-Sugimoto 
model have been computed in [13] (see also [14,15]), by performing a "small-size" {i.e. small 
7) expansion in which analytical results can be obtained. This expansion is not trustable 
in our case because, as explained above, the size of the 5D skyrmions is large and 7 ~ 1. 

The paper is organized as follows. In sect. 2, after a brief review of the model and of 
the static 5D skyrmion solution we identify the zero-mode fluctuations which are relevant to 
describe static properties and we discuss the corresponding collective coordinates classical 
Lagrangian. A suitable ansatz is described which permits to rewrite in a 2D form the 
4D equations which define the zero-modes. Sect. 3 is devoted to the collective coordinate 
quantization and to the calculation of the form factors, this discussion is basically the same as 
in the 4D Skyrme model [6,7], though adapted to the present case. Sect. 4 contains a detailed 
presentation of our results. After the comparison with simple expectations about laige-Nc 
QCD we discuss the divergences due to the chiral limit and we check that the Goldberger- 
Treiman relation holds in our model. Finally, we compare our results with experimental 
data and find a level of agreement better than 30% for all the observables, with the notable 
exception of the axial coupling gj^ for which we find qa = 0.70 versus an experimental value 
gA = 1.25]^ Most of the technical details are presented in the appendices. In appendix A 
the 2D equations of motions and boundary conditions are derived while appendix B gives 
some detail on the numerical techniques we employed to obtain the solution. 

In spite of the failure in the axial coupling, the level of accuracy of our results is consistent 
with the expected size of the 1/Nc corrections or, which is the same, the expected size of 
next-to-leading contributions in our model. It is not unreasonable that anomalously large 

^The same calculation has been performed in [16] with a different (and erroneous, in our understanding) 
definition of the chiral currents. 

^An erroneous value of gA was reported in [9]. The error was due to a subtlety, which we will discuss 
in the following, in taking the zero momentum limit of the axial form factor, combined with a more trivial 
mistake. 
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numerical factors could change into 80% the naively expected 30% correction to qa- Such 
large 1/Nc corrections arise for example if one follows, instead of the approach we consider, 
the quantization procedure of the collective coordinates proposed in [17]. This "alternative" 
quantization is equivalent to the standard one at the leading order in l/N^. but it also 
contains large l/N^ corrections. We will discuss in sect. 4 how these corrections change our 
predictions. In the case of qa we find, remarkably, the much better result qa = 1-17 while 
the level of the agreement of the other observables is unaffected. Even without applying 
this correction, the results which we obtain are significantly more accurate than those of the 
original Skyrme model [6] (in which, we remark, qa is also small, qa = 0.65), but not as 
good as those of more refined skyrmion models (which, of course, also have more parameters) 
like the ones reviewed in [7]. It seems, as we will discuss in sect. 4, that the inclusion of 
the explicit breaking of the chiral symmetry {i.e. of the pion mass m^r) will improve the 
agreement of several observables and that qa may display an enhanced sensitivity to m^. It 
is certainly worth exploring this direction. 



2 Skyrmions in 5D 
2.1 The Model 

We will consider the same model as in [9], i.e. a U{2)l x U{2)r gauge theory in five 
dimensions with metric ds^ = a(z)^ {rj^ydx^dx'^ — dz"^), where we denoted as x'^ the usual 4 
coordinates with mostly minus metric and with z, which runs in the interval [2;uv,-2ir], the 
extra dimension. We choose our metric to be AdSs and therefore the warp factor a{z) to be 



with z\jY — to be taken at the end of the calculations. In this limit, zj-r coincides with the 
conformal length L = dz = — zyy. It should be kept in mind that, since gravity is 
non-dynamical in our model, the choice of the warp factor a{z) is arbitrary. It is commonly 
believed, however, that different "reasonable" choices of a{z) would not affect in a significant 
way the predictions for IR observables, like those we will compute in this paperj^ Choosing 
AdSs -or at least a geometry with an "AdSs-like" boundary- is crucial, on the contrary, if 

■^We thank the referee of Nucl. Phys. A for suggesting this possibihty to us. 

^This behef is supported by Ref. [18], in which IR predictions for flat and AdSs spaces were compared, 
and by Ref. [3], in which departures from AdSs were considered. Moreover, the Sakai-Sugimoto model is 
equivalent, for what calculations in the meson sector are concerned, to a 5D model of the kind we are 
considering with non-AdSs warped metric. The predictions of this model are very similar to those of 
AdS/QCD, again suggesting that the choice of the metric is not so relevant. 
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one wants to match UV correlation functions with those computed in QCD by perturbation 
theory [1-4,18]. For this reason, in the hterature the choice in eq. ([T]) is commonly adopted. 

We will denote the U{2)i and f/(2)/j gauge connections respectively by and Rm, 
where M = {jj, 5}, and parametrize them as = L1jaa/2 + Lm^/'^ and Ra/ = R%jaa/2 + 
Rm^/2 in terms of the Pauli matrices (Ja and the identity 1. Chiral symmetry is broken at 
the z = zjji boundary (IR-boundary) by the following conditions: 

(L^ - R/^) \z=zm = ) + R-ms) |2=zir = , (2) 

where the 5D field strength is defined as Lmn = c^mLtv — c^AfLM — ^[Lm, Ltv], and analogously 
for TLmn- On the other boundary, the UV one, we impose Dirichlet conditions: 



The 5D action S = Sg + Scs consists of a standard gauge kinetic part 
S9 = -l d'xj^^'dzaiz) ^ jxr [LmnL"""] + ^LmnL^'" + {L ^ , (4) 

and of a Chern-Simons part 

Scs = ^,j rf^x|le*^^^^«LMTr[L^o^PQ] + ^e^^^^^^L^L^o^Pg - {L ^ R]^. 

(5) 

The Scs is needed to reproduce the QCD anomalies and its coefficient is fixed to be propor- 
tional to the number of colors Nc- 

In order to compare our 5D model with the real world, and in particular to compute the 
form factors as we will do in this paper, we need to identify the chiral currents to which 
the electroweak bosons are coupled. These operators, which would be given in QCD by the 
quark bilinears 31,^l{r) = Ql(R)1^(^'' I'^Ql{R), 3tM,L{R) = Ql{r)Y^I2Ql{R) correspond in our 
model to [9] 

Jl^ = M^{a{z)Ll,) 1,=,^, , J,^ = a^M^{a{z)L^,) , (6) 

and analogously for R. 

It is important to remark that this model, as discussed in [9], is a valid effective field 
theory with an NDA cut-off A5 which is bigger than the scale of the lightest resonances and 
can be sent to infinity for M5 — 00. This allows us to include only the lowest dimensional 
operators in the action 



ifS ) since the others, which would surely arise in a UV completion of 
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the model, are expected to give a subleading contribution. At the leading order this model 
is extremely predictive: its only 3 parameters are M5, L and a. The 5D interaction scale 
M5 can be traded for the adimensional parameter 

(7) 

which controls the size of the skyrmion; p ~ 7^/^. We want to interpret the 5D weak coupling 
expansion as the l/N^ expansion, therefore we will take the interaction scale to scale like iVc, 
i.e. M5 oc Nc so that 7, p oc N'^. 

2.2 The Static Soliton Solution 

Our model admits topologically non-trivial static solutions of the classical equations of mo- 
tions (EOM). These are identified with the baryons and therefore the topological charge 

- I fl^x I 

327r2 

is identified with the baryon number. The indices /i, z>, . . . label, throughout the paper, the 
4 spatial coordinates, but they are raised with Euclidean metric. 



Regular static solutions with B = 1 have been found in [9]. The non- vanishing compo- 
nents of the R fields can be written in terms of 2D fields as 

Rl{^,z) = A2{r,z)x\ (9) 
% I ^ s{r, z) 
r 

where = = x^/r, is the antisymmetric tensor with e^-^'^^ = 1 and the 

"doublet" tensors A*^-*^'^) are 



^{x),ab 



(10) 



Due to parity invariance {L ^ i?, x — x} we restrict, in both the static and non-static 
case which we will consider in the next section, to configurations for which Lj(x, = 
—Ri{—x,z,t), L5fl{x,z,t) = R5fl{—x,z,t) and analogously for L, R. Eq. ^ therefore 
defines the static solution completely. 



It is important to remark that the static solution in eq. (|9j) is "cylindrically" symmetric, 
meaning that it is invariant under the simultaneous action of 3D space rotations XqCt" 
r'^XaO'°'r, with r G SU{2), and vector SU{2) global transformations L,R^ r{L,R)r''. An 
equivalent way to state this is that a 3D rotation with r acts on the solution (|9| exactly as 
an SU{2) vector one in the opposite direction {i.e. with r"!") would do. 
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2.3 Zero-Mode Fluctuations 



Let us now consider time-dependent infinitesimal deformations of the static solutions. Among 
these, the zero-mode {i.e. zero frequency) fluctuations are particularly important as they 
will describe single-baryon states. Zero-modes can be defined as directions in the field space 
in which uniform and slow motion is permitted by the classical dynamics and they are 
associated with the global symmetries of the problem, which are in our case U{2)v and 
3-space rotations plus 3-space translations. The latter would describe baryons moving with 
uniform velocity and therefore can be ignored in the computation of static properties like 
the form factors. Of course, the global U{l)v acts trivially on all our fields and the global 
SU{2)v has the same effect as 3-space rotations on the static solution (joj) because of the 
cylindrical symmetry. The space of static solutions which are of interest for us is therefore 
parametrized by 3 real coordinates -denoted as collective coordinates- which define an SU{2) 
matrix U. 

To construct zero-modes fluctuations we consider collective coordinates with general time 
dependence, i.e. we perform a global SU{2)y transformation on the static solution 



i?/i(x, z; U) 



URf,{^, z) 



Ro{yi,z]U) = -Ro(x, 



'111 



but we allow U = U{t) to depend on time. It is only for constant U that eq. ( 11 ) is a solution 



of the time-dependent EOM. For infinitesimal but non-zero rotational velocity 

K = kaay2 = -iUUU/dt , 



eq. (11) becomes an infinitesimal deformation of the static solution. Along the zero-mode 



direction uniform and slow motion is classically allowed, for this reason our fluctuations 
should fulfill the time-dependent EOM at linear order in K provided that dK/dt = 0. 



From the action (4][5) the following EOM are derived 



Do (a(^)/?y + 



■jaL 
4 

ado {a{z)R%) + 



.vuipa 



RoujR 







MLUpa 



TriRo^R 



Do {a{z)R^f') - a{z)DoRo^ - 



■yaL 



pa J 



pvpa 



. R'ui'iiR^ 



RuoRpa + RuoR, 



ado [a{z)R^n - aa{z)doRQ^ - -fLe 



fivpa 



:i2) 



Tr (RooRpa) + ^RuoRpa 
We only need to specify the EOM for one chirality since we are considering, as explained in 



the previous section, a parity invariant ansatz. We would like to find solutions of eq. (12) 
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for which i?^ and Ro are of the form (11); it is easy to see that the time-dependence of 



[/ in eq. (11) acts as a source for the components Rq and Rn, which therefore cannot be 



put to zero as in the static case. Notice that the same happens in the case of the 4D 
skyrmion [7], in which the temporal and spatial components of the p and u mesons are 



turned on in the rotating skyrmion solution. Also, it can be shown that eq. (12) can be 
solved, to linear order in K and for dK/dt = 0, by the ansatz in Eq. (11) if the fields Ro 



and RfX are chosen to be linear in K. Even though K must be constant for the EOM to 
be solved, it should be clear that this does not imply any constraint on the allowed form of 



the collective coordinate matrix U{t) in eq. (11), which can have an arbitrary dependence 
on time. What we actually want to do here is to find an appropriate functional dependence 
of the fields on U (t) such that the time- dependent EOM would be solved if and only if the 
rotational velocity K = —iWdU /dt was constant. 



In order to solve the time-dependent equations (12) we will consider a 2D ansatz obtained 



by a generalization of the cylindrical symmetry of the static case. The ansatz for i?^ and Rq 
is specified by eq. (11) in which the static fields are given by eq. (|9]). Due to the cylindrical 



symmetry of the static solution the fields in eq. (|11|) are invariant under 3D space rotations 



> r''XaCr°'r combined with vector SU (2) global transformations L,R —>■ r (L, R) r'^ if 
U also transforms as U ^ r'^Ur. We are therefore led to consider a generalized cylindrical 
symmetry under which ka also rotates as the space coordinates do. Compatibly with this 
symmetry and with the fact that Rq and i?^ must be linear in K we write the ansatz as 



Ro{x,z;U) = U Ro{-K,z;K)U^ + iUd^U^ 



(13) 



where 



rI{^,z-K) = 
aRii^, z; K) 



X{x)[.r,z 
p{r, z) 



)ki,A^-)''''^ + v{r,z){k-x)x^ 

[k' - {k ■ x)x^) + 5i(r, z){k ■ x)x^ + Q(r, z)e^^''hx. 



(14) 



ai?5(x, z] K) = B2{r, z){k ■ x) 



It should be noted that the term iUdoW = UKU'^ in eq. (13) is purely conventional 



as it could have been reabsorbed in the definition of Rq. Nevertheless this choice makes 
manifest that our ansatz (TT|l3) can be obtained from the "barred" fields in eq. ( 9p4 ), 
which only depend on U through K., by performing a time-dependent SU{2) vector gauge 
transformation with parameter U{t). This is useful because the action, including the CS 
term, is invariant under this transformation. We can therefore obtain the 2D EOM for our 
ansatz fields by plugging the barred fields, instead of the original ones, into the 5D EOM. 
It is important to stress that the ansatz with barred fields is not truly gauge equivalent to 
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the original one because the transformation U (t) does not reduce to the identity at the UV 
boundary, implying that the UV condition (|3| is not invariant. Our true ansatz is therefore 
provided by eq. ( ITp3 ) and the use of the barred field as we will do in the following is just 



a useful trick. 

At this point it is straightforward to find the zero-mode solution. The EOM for the 2D 



fields can be obtained by plugging the ansatz in eq. (12), while the conditions at the IR and 
UV boundaries are derived from eq. (|2| and (|3]), respectively. The boundary conditions at 
r = are obtained by imposing the regularity of the ansatz, while those for r — > oo come 
from requiring the energy of the solution to be finite and the topological charge B in eq. ^ 
to be equal to 1. More details are presented in appendix A, where the 2D EOM and the 
boundary conditions are derived. Once the 2D equations have been found, however, it is not 
yet trivial to solve them numerically, the procedure we followed is described in appendix B. 



The reader not interested in detail can however simply accept that a solution of eq. (12) 
exists and is given by our ansatz for some particular functional form of the 2D fields which 
we are able to determine numerically. In the rest of the paper the 2D fields will always 
denote this numerical solution of the 2D equations. 



2.4 The Lagrangian of Collective Coordinates 



The collective coordinate matrix U(t) will be associated with static baryons. The classical 
dynamics of the collective coordinates is obtained by plugging eq. ( TTp3 ) in the 5D action. 
One finds ^[f/] = J dtL where 

L = -M+^A:,fcr (15) 
The mass M and the moment of inertia A are given respectively by 



M = SttM. I dr dz i a(z) 



•yL s 

( 

2 r 



4 



d-^{-i(f)*Dp(j) + h.c.) + A 



flU 



and 



^ POO PZl^ 

A = 167rM5- dr dz { aiz) 



-J 2UV 



- {D-^pf - r^ {d-^Qf - 2Q2 _ 



(Dfix) + y (dfiv) + {x(x)X(x) + v^) (l + (f)(x)(f)(x)) - ^v(f)(^^)X{x) 



2e^^D-^pXix) {D, 
Qe^^E^p (0(.)0(x) - l)+rQef''Af,^^ + 



2e-^'dp{rQ) X(x)e'^''^ (I^p0)(,) 



(16) 



(17) 
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The notations used in the equations above are defined in appendix A; the covariant derivative 
symbols, in particular, are associated with two Abelian residual gauge symmetries which our 
2D ansatz has. Here we simply want to show that M and A could be easily computed, at a 
given point of the parameter space, once the numerical solution for the 2D fields is known, 
by performing a numerical 2D integral. 

Let us give some more detail on this theory. For now we proceed at the classical level 
and we will discuss the quantization in the next section. Our Lagrangian can be rewritten 
as 



L = -M + ATr 



where we have parametrized the collective coordinates matrix U a.s U = UqH + iuia'^, with 



^ . = 1. The Lagrangian ( 18 ) is the one of the classical spherical rigid rotor. The variables 
{uo,Ui} are restricted to the unitary sphere S^, which is conveniently parametrized by the 
coordinates g° = {x, 0i,02} -which run in the x G [—1,1], 0i G [0,27r) and 02 S [0,27r) 
domains- as 



u,+tu2 = zi = ^^^e"^\ uo + tU3 = Z2 = ^^^e'^\ (19) 

where we also introduced the two complex coordinates Zi^2- We can now rewrite the La- 
grangian as 

L = -M + 2Xg^pq''f , (20) 
where g is the metric of which reads in our coordinates 

ds^ = gapdq^dq^ = \-^-^dx^ + ^—^ d^i^ + • (21) 

4 1 — 2 2 

The conjugate momenta are Pa = dL/dcf = AXg^fsq^ and therefore the classical Hamiltonian 
is 

H, = M+ l^p^g-P{q)pp. (22) 

It should be noted that the points U and —[/in what we denoted as the space of collective 
coordinates actually describe the same field configuration (see eq. ( Ilp3 )). The SU{2) = 5; 
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manifold we are considering is actually the universal covering of the collective coordinate 
space which is given by 83/ Z2. This will be relevant when we will discuss the quantization. 
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3 Static properties of Nucleons 



3.1 Skyrmions Quantization 



We should now quantize the classical theory described above, by replacing as usual the 
classical momenta pa with the differential operator —id/dq" acting on the wave functions 
/(g). Given that the metric depends on q, however, there is an ambiguity in how to extract 



the quantum Hamiltonian Hg from the classical one in eq. (22). This ambiguity is resolved 



by requiring the quantum theory to have the same symmetries that the classical one had. 
At the classical level, we have an S'0(4) ~ SU(2) x SU{2) symmetry under U U ■ r'^ and 
U ^ g ■ U with r,g E SU{2). These correspond, respectively, to rotations in space and to 
isospin (i.e. global vector) transformations, as one can see from the ansatz in eq.s ( Tlp3 ). 
This is because K is invariant under left multiplication by g, and that the ansatz is left 
unchanged by performing a rotation XaC"" — > r'^XaC'^r and simultaneously sending U —>■ U -r. 
The spin and isospin operators must be given, in the quantum theory, by the generators of 
these transformations on the space of wave functions f{q) which are defined by 



After a straightforward calculation one finds 



(23) 



5+ 

s- 

/+ 



1 

1 

i 



iVl- x'^dx + l:\/T——dcf>i - — -dd, 



2 V 1 - X 



2 V 1 +a; 



-i{<f)i+4>2) 



1 



2 V 1 -X 



2 \ 1 + x 



(24) 



2 V 1 -X 



2 V 1 + x 



-i{(t>l-(t>2) 



«Vl - x^dx - -\j d^-, - 

X 



2yi-x''' 2Vl + 



where the raising/lowering combinations are S*^ = {S^ ± iS'^) / \/2. 



The operators in eq. (24) should obey the Hermiticity conditions (5^)^ = 5'^, (5*+)^ = S 



and analogously for the isospin. In order for the Hermiticity conditions to hold we choose 
the scalar product to be 

{A\B) ^ [ d\^fA\q)fB{q), (25) 
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where ^Jg = 1/4 in our parametrization of S^. The reason why this choice of the scalar 
product gives the correct Hermiticity conditions is that 5" and /" (where a = 1,2,3) can 
be written as X'^da with Kilhng vectors of the appropriate isometries. The Kilhng 
equation V^X^ + V/jXq, = ensures the generators to be Hermitian with respect to the 
scalar product (25). 



Knowing that the scalar product must be given by eq. (25) greatly helps in guessing 
what the quantum Hamiltonian, which has to be Hermitian, should be. We can multiply 
and divide by y/g the kinetic term of He and move one y/g factor to the left of pa- Then we 
apply the quantization rules and find 

H, = M- ^^5a i^gg^'^dp) = M- ^V<,V" , (26) 

which is clearly Hermitian. We can immediately show that Hq commutes with spin and 
isospin, so that the quantum theory is really symmetric as required: a straightforward cal- 
culation gives indeed 

Hq = M + = M + . (27) 



It would not be difficult to solve the eigenvalue problem for the Hamiltonian (26), but in 
order to find the nucleon wave functions it is enough to note that the versor of n-dimensional 
Euclidean space provides the n representation of the SO{n) isometry group. In our case, 
n = 4 = (2,2), which is exactly the spin/isospin representation in which nucleons live. It 
is immediately seen that zi, as defined in eq. (19), has = = 1/2. Acting with the 
lowering operators we easily find the wave functions 



1 i 

\P T) = -zi, \n T) = -Z2, 

\P i) = Z2, \n i) = Zi 

71 71 



{2i 



which are of course normalized with the scalar product (25). The mass of the nucleons is 
therefore E = M + 3/(8A). 

Notice that the nucleon wave functions are odd under U —U, meaning that they are 
double- valued on the genuine collective coordinate space 83/22. This corresponds, following 
[19], to quantize the skyrmion as a fermion and explains how we could get spin-1/2 states 
after a seemingly bosonic quantization without violating spin-statistic. 

Let us now summarize some useful identities which will be used in our calculation. First 
of all, it is not hard to check that, after the quantization is performed the rotational velocity 



'^The last equality holds because Hq is supposed to be acting on the wave functions, which are scalar 
functions. 



11 



becomes 



and analogously 



r = -iTi 



i Tr 



A 



(29) 
(30) 



Other identities which we will use in our calculations are 



(Tr [f/(T^f/V"] = --s^r) 



(Tr [Ua^Xb{k ■ x)U^a^] = -^/") 



(31) 



where the VEV symbols (...) mean that those are not operatorial identities, but they only 
hold when the operators act on the subspace of nucleon states. Notice that the second 



equation in (31) is implied by the first one if one also uses the commutation relation (23) 



eq. (30) and the fact that, on nucleon states, ({5"", S^} = 6'^^/2). 



3.2 The Nucleon Form Factors 



The nucleon form factors parametrize the matrix element of the currents on two nucleon 
states. For the isoscalar and isovector currents we have 



{Nfip')\4{o)mp)) = ujip') 



(iV/(p')l^r(o)|iv.(p)) = ^/(p') 



2Mn 



2M, 



N 



(32) 



where the currents are defined as Jy = J'^ + J£ and Js = 1/3 (^Jr + J^j in terms of the 
chiral ones. In the equation above q = p' — p is the 4-momentum transfer, Ni and Nf are 
the initial and final nucleon states and Ui{p), Uf{p') their wave functions, J" = a"/2 is the 
isospin generators and a^'^ = i/2['j^,j'^]. For the axial current Ja = Jr^ Jl have 



{Nf{p')\J'X^mN,{p))=Uf{p')GA{q') 



2M, 



7m 



N 



-q' 



ri^ufip) . 



(33) 



Exact axial and isospin symmetries, which hold in our model, have been assumed in the 
definitions above. 

In our non-relativistic model the current correlators will be computed in the Breit frame 
in which the initial nucleon has 3-momentum —q/2 and the final +q/2 (i.e. p'^ = {E, —q/2) 
and p''^ = {E,q/2), and = — g^, with E = 



'Mli + q'^/A). Notice that the textbook 
definitions in eq.s ( 32||33 ) involve nucleon states which are normalized with \/2E\ in order 
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to match with our non-relativistic normahzation we have to divide all correlators by 2Mn. 
The vector currents become 

(iV;(g72)|J°(0)|iV,(-g72)) = GUf)x\x^, 

{Nf{q/2)\JhmN^{-m) = ^ ^^^2(5^ X , 
{Nf{q/2)\J'viO)m-q/2)) = GUf)x\{2nX^, 

{Nfiq/2)\Jmm-q/2)) = ^^f^x\2{S x q)^ {2n x^ , (34) 
where we defined 

Gn-,') = Ff'"(g^) + , G^/i-q^) = Ff'-(g^) + Ff'-(g^) , (35) 

and used the definition (5 x g)* = e''^^S^q^. The nucleon spin/isospin vectors of state Xij 
are normalized to x^X = 1- Fo^' the axial current we find 

U)|7V,(-g72)) = x}-^ 
{Nf{q|2)\JTmN^{-q|2)) = Q (36) 
where St = S — q S ■ q is the transverse component of the spin operator. 



{Nj{ql2)\JTmN^{-ql2)) = x\J;r^A{f)2S!,^x^ 



It is straightforward to compute the matrix elements of the currents in position space on 
static nucleon states. Plugging the ansatz ( [9|ll[14[13[ ) in the definition of the currents ^ 
and performing the quantization one obtains quantum mechanical operators acting on the 
nucleons. The matrix elements are easily computed using the results of sect. 3.1. We finally 
obtain the form factors by taking the Fourier transform and comparing with eq.s ( [34p6| . 
We have 

Nc f 

Ge = J drrjo{qr){a{z)d,s)^y 

Gl = ^ I drr'Uqr) [a{z) (d^v - 2{D^x)i2) 



uv 



E Sna'jL 



jiiqr) ( \ 
a{z) — — \ {p4)^^■^ -rA.rJ -a{z) {D,4>)^^^jo{qr) 



(37) 

uv 



JiiSF) 
qr 

where j„ are spherical Bessel functions which arise because of the Fourier transform. 

^It is quite intuitive that the form factors can be computed in this way. Given that sohtons are infinitely 
heavy at small coupling, in the Breit frame they are almost static during the process of scattering with the 
current. To check this, however, we should perform the quantization of the collective coordinates associated 
with the center-of-mass motion, as it was done in [20] for the original 4D Skyrme model. 
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4 Results 



In this section we will present our results. After discussing some qualitative features, such 
as the large- Ai'c scaling of the form factors and the divergences of the isovector radii due 
to exact chiral symmetry, we extrapolate to the physically relevant case of A^c = 3 and 
perform a quantitative comparison with the experimental data. Consistently with our work- 
ing hypothesis that the 5D model really describes large-iVc QCD we find a 30% relative 
discrepancy. 



Large- Ac Scaling 

Let us take all the three parameters a, 7 and L of our 5D model to scale like A^ for large- A^. 
Eq. ([7]) therefore implies that the coupling M5 grows like A^ and the semiclassical expansion 
in 5D coincides with the 1/Nc expansion on the 4D side. Notice that these scaling of the 
parameters are uniquely dictated by what we know to be the large-A"c scaling of meson 
couplings and masses. In the baryon sector, the solitonic solution is independent of N^. given 
that M5 factorizes out of the action and does not appear in the EOM. The classical mass M 
and the moment of inertia A therefore scale like and the scaling of the form factors can 



be easily read from eq. (37) 



In large- A"c QCD the scaling of several baryon observables is known [21]. The mass grows 
with Nc as in our model, but this is expected to be a common feature of any soliton model. 
The matrix element of currents on normalized nucleon states should be of the form N^F{q'^) 
with p = 1 even though cancellations, i.e. p < 1, are not excluded. All the radii should 
therefore scale like A^ and this is what we find in our model. We also find the "naive" -i.e. 
with p = 1- overall scaling for the electric scalar (Gf), magnetic vector {G\.f) and axial (Ga) 



form factor; notice that, due to the definition in eq. (34), the magnetic form factors scale 
with one more power of than what the corresponding current matrix element does. We 
however find two cancellations: due to the 1/A factor the electric vector and the magnetic 
scalar Gf,j scale like and N^, respectively. This corresponds to a "p = 0" scaling of the 
associated currents. 

The reason for the cancellation in G^ is very simple to understand. Remembering that the 
temporal component of the current at zero momentum gives the conserved charge and looking 



at the definitions (34) one immediately obtains two consistency conditions: G^{0) = N^/G, 
because in the nucleons there are Nc quarks which have U{l)v charge 1/6 each in our 
conventions, and G^{0) = 1/2, because nucleons are in the 1/2 representation of isospin. It 
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is not difficult to see that these consistency conditions are respected by our model as they 
are implied by the EOM, and they are fulfilled to great accuracy (0.1%) by the numerical 
solution. The above discussion implies, in particular, that while the electric scalar form 
factor G% has the naive scaling, the electric vector does not. 

We are not able to prove that the cancellation in Gfj actually takes place in large-A^^c 
QCD, but we can check that it occurs in the naive quark model, or better in its generalization 
for arbitrary odd = 2 k+1 [22]. In this non-relativistic model the Nucleon wave function is 
made of 2A;+1 quark states, 2k of which are collected into k bilinear spin/isospin singlets while 
the last one has free indices which give to the Nucleon its spin/isospin quantum numbers. 
Of course, the wave function is symmetrized in fiavor and spin given that the color indices 
are contracted with the antisymmetric tensor and the spatial wave function is assumed to 
be symmetric. The current operator is the sum of the currents for the 2k + 1 quarks, each of 



which will assume by symmetry the same form as in eq. (34). If Si^2 and Ii^2 represent the 
spin and isospin operators on the quarks qi^2 the operators 5*1 + 5*2 and Ji + I2 will vanish 
on the singlet combination of the two quarks, but Sili + S2I2 will not. The k singlets will 
therefore only contribute to , G\j and Ga, which will have the naive scaling, while for 
the others we find cancellations. 

A detailed calculation can be found in [23] where, among other things, the proton and 
neutron magnetic moments and the axial coupling are computed in the naive quark model. 
The magnetic moments are related to the form factor at zero momentum as ^v/f^N = G'm(O) 
and Hs/lJ'N = CfflO) where /itv = l/(2Mjv) is the nuclear magneton and 2 = fip — fin, 
2 Us = fip + fin- In accordance with the previous discussion, the results in the naive quark 
model are 2/is = fiu + fid and 2/iy = 2/c/3(/i„ — fid), where fiu,d are the quark magnetic 
moments, while for the axial coupling one finds qa = Ga{0) = 2k/3 + 1 = N^/S + 2/3 which 
scales like as expected. Notice that the 1/Nc corrections to the axial coupling are quite 
big in the naive quark model: for Nc = 3 the leading term contributes as 1 while the "true" 
results is 5/3, which is 67% bigger. We have of course no reason to believe that such big 
corrections should persist in the true large- A*";, expansion of QCD, this trivial remark simply 
suggests that "large" l/A'c corrections to the form factors are not excluded. 



Divergences in the Chiral Limit 

It is well known that in QCD the isovector electric (r|; y) and magnetic (r|j y) radii which are 
proportional, respectively, to the derivative of G^ and G\j at zero momentum, diverge 
in the chiral limit [24]. If the nucleons are effectively described by an isospin doublet of 
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point-like spinors added to the xpt Lagrangian this effect comes from pion loops which are 
IR divergent in the massless pion limit —>■ [25]. It is not obvious, however, that the 
divergences should survive in the large- A^c limit, i.e. that they should already appear in the 
leading term of the perturbative 1/Nc expansion. This is so because describing baryons as 
weakly coupled particles, which is a reasonable approximation in real-world Nc = 3 QCD, 

3 /2 

is not possible at large- iVc given that their coupling with pions grows like Nc . It was 
noticed in [26] that, in a model which only contains the nucleon doublet, the pion-nucleon 
scattering amplitude grows like Nc violating unitarity and also contradicting the usual large- 
Nc counting rules. The theory is therefore inconsistent and the full infinite tower of large- iVc 
baryons must be added. Moreover, it is very easy to see that the results of [24, 25] are not 
compatible with large Nc- the one-loop corrections to the radii, which are of course finite for 
finite m^, have the wrong scaling and grow like Nc- Given that we cannot apply the results 
of [24, 25] , we cannot conclude that the radii must diverge in our model at the leading order 
in the semiclassical expansion of the soliton, but there is of course no problem if they do. 
We must however check that all the other radii are finite, and this is what we will do in the 
following. What we will find is the same as in the 4D Skyrme model [6]: all radii and form 
factors are finite but the electric and magnetic isovector ones. 



In our model, as in the Skyrme model, divergences in the integrals of eq. (37) which 



define the form factors are due, as in QCD, to the massless pions. If all the fields were 
massive, indeed, any solution to the EOM would fall down exponentially at large r while 
in the present case power-like behaviors can appear. These power-like terms in the large-r 
expansion of the solution can be derived analytically by performing a Taylor expansion of the 
fields around infinity (1/r = 0), substituting into the EOM and solving order by order in 1/r. 
The exponentially suppressed part of the solution will never contribute to the expansion. In 
the gauge in which (the form factors are, of course, gauge invariant) the topological twist is 
at the origin r = and the solution is trivial for r 

2z{z - 1) 



A2 



z(l -z)^^ z{z^-2z^ + l) ^ 



oo the first few terms are 
ziz-1) 



■1 + 



z^{3-2z 



2r^ 



2r^ 
Az^+ 



4^8 



Xi 

X2 
V = 



1 



'z\2z-3) 



-1 + 
zHz^ 



z^iz^ - 



3f 



12r^ 
-Az^ + 



4^8 



{31 



where f3 is an unknown parameter which depends on the entire solution and can only be 



^In the equations which follow we put L ~ 1 for simplicity. 
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determined numerically. We checked that the large-r behavior of our numerical solution is 



very well approximated by eq. (38). Substituting these expressions into the definitions of 



the form factors (37) one gets 



Gl oc (3' 



dr^joiqr) + ... 



Gl ^ (3^ I dr-Uqr) + ... 



dr\'^ + ... 



(39) 



Gl, cx 



dr 



J.2 qj. 



+ 



All the form factors are finite for any q, including g = 0. The electric and magnetic radii, 
however, are defined as 



6 



dGE,M{q^] 



GE,M{q'^ 



0) dq- 



(40) 



and taking a q^ derivative of eq.s (39 ) makes one more power of appear in the integral. It is 



easy to see that the scalar radii are finite, while the vector ones are divergent as anticipated. 
We will now discuss the axial coupling and the axial radius and show that both are finite. 



The expression in eq. (37) for the axial form factor Ga presents some subtleties for 
vanishing q^. Given the asymptotic expansion of the solution in eq. (38) the axial coupling 
integral behaves for large r like 



Ga (X dr 



:/3 - -J' 



ji{qr) 



qr 



r 



5 

7r^ 



) Mqr) + ... 



(41) 



where the leading 1/r terms (the ones which are linear in (3) can be obtained from eq. (38) 



while for the others one needs higher order terms which are not reported in eq. (38). The 



integral in eq. (41) is finite integral for any q ^ 0. For g — 0, however, the integral is not 



uniformly convergent and one cannot exchange the limit with the integration. The leading 



1/r term in eq. (41) is indeed given by J(g) = P dr (1/r) (3ji(gr)/(gr) — jo{qr)) which 



is independent of q and equal to /?/3. Given that the argument of the integral vanishes for 
q ^ exchanging the limit and integral operations would give the wrong result J(0) = 0. 
To restore uniform convergence and obtain an analytic formula for qa one can subtract the 



I{q) term from the expression in eq. (37) for Ga and replace it with /?/3. Rewriting the axial 
form factor in this way is also useful to establish that the axial radius, which seems divergent 



if looking at eq. (41), is on the contrary finite. The I{q) term, indeed, does not contribute 



to the derivative and the ones which are left in eq. (41) give a finite contribution. 
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We have found, in summary, that all the form factors and radii are finite but the isovector 
ones. Given that the divergences are related with the large-distance behaviors of the fields, 
and that our model reduces to the Skyrme model in the IR, this result is not surprising. A 
different result has been found, however, in Ref. [13], where the nucleon form factors have 
been computed in the Sakai-Sugimoto model. In that case all the radii are finite. The Sakai- 
Sugimoto baryons correspond, as explained in the Introduction, to the small-size limit of the 
5D skyrmions we are considering, and we should recover the results of [13] if we perform a 
small-7 expansion which correspond to the 1/A expansion considered in [13]. As 7 decreases 
our soliton becomes more and more localized around {r = 0,z = ^ir) and at any large 
but fixed value of r the deviations from the pure-gauge configuration become smaller and 



smaller. The small-7 expansion of the asymptotic solution (38) therefore coincides with the 



small-/? expansion. By looking at eq. (39) we see that the power-like terms in the isoscalar 
and isovector form factor densities appear at high orders in (3 and this explains why these 



densities were found to be exponentially damped in [13]. For the axial form factor, as eq. (41 ) 
shows, power-like terms are present at the linear order in /3. The same term has been found 
in [13] but it does not lead to any divergence as explained in the previous paragraph. 

A possible physical explanation of the finiteness of the radii in the Sakai-Sugimoto model 
is that the 5D soliton effectively reduces to a 5D particle in the limit of small size, a possibility 
discussed in [28]. For a 5D particle no divergences appear in the radii at the leading order in 
the semiclassical expansion {i.e. at tree level) and the divergences should arise, in analogy 
with the case of a 4D particle, at loop level. Following the analogy, however, one could 
expect the divergent loop corrections (or better the enhanced loop corrections for small but 
finite pion mass) to have, as it happens for the 4D particle, the wrong large-A^^c scaling. 
By the same reasoning one could expect unitarity violation in the pion-nucleon scattering 
amplitude at tree-level. 



Pion Form Factor and Goldberger— Treiman relation 

It is of some interest to define and compute the pion-nucleon form factor which parametrizes 
the matrix element on Nucleon states of the pion field. In the Breit frame (for normalized 
nucleon states) it is 

(Ar;(g72)|7r'^(0)|Ar,(-g72)) = __^G;v^^(g-«);^t (2^0g,(2r)x. , (42) 

where n°-{x) is the normalized and "canonical" pion field operator. The field is canonical 
in the sense that its quadratic effective Lagrangian only contains the canonical kinetic term 
£2 = l/2(97ra)^, or equivalently that its propagator is the canonical one, without a non-trivial 
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form factor. With this definition, Gnntt is the vertex form factor of the meson-exchange 
model for nucleon-nucleon interactions [29] and corresponds to an interaction 

JO-NNn = I (G^iv.(n)vrJiV7'^75(2/")iV . (43) 

On-shell, the form factor reduces to the pion-nucleon couphng constant, Gnntt{S^) = dNNn, 
whose experimental value is gNNn = 13.5±0.1. 

The pion field which matches the requirements above is given by the zero-mode of the 
KK decomposition. In the unitary gauge dz{a{z)A^) = 0, where Am = {Lm — Rm)/2, and 
for AdSs space a{z) = L/z one has 

with = 2M5/ J dz/a{z) = AM^/L. |^ Gauge-transforming back to the gauge in which our 
numerical solution is provided and using the ansatz in eq.s dgfll] ) we find the pion field 



vr 



2 



/ rf2A2(r,2)J*Tr [[/(Tfcf/V'^] . (45) 



Taking the matrix element of the above expression and comparing with eq. (42) one obtains 
Gjviv.(g') = -yM^F^gy drn{qr) j dzr^ A^ir.z) . (46) 



Using eq. (38) it is easy to understand that the g — >■ limit of Gnn-k is completely 
determined by the large-r behavior of the field y42, and in particular by the leading 
term. Due to the q factor, indeed, only the divergent part of the integral contributes. We 
then find 

gNNn = — :r- 13 . 47 

We used the formula above to check numerically that the Goldberger-Treiman relation 
F-kQ-kNN = Mj^qa holds in our model, we find that it is verified to 0.01% on our numer- 
ical solution. We can also demostrate the Goldberger-Treiman relation by using eq. (33,34) 
of [6] which show that also qa is determined by the asymptotic behavior of the axial current. 
We indeed obtain 

9. = -|^/J. (48) 



^Nucleon scattering, in our model, is a soliton scattering process and we have no reason to believe that it 
can be described by meson-exchange, i.e. that contact terms are suppressed. Therefore, we will not attempt 
any comparison of our form factor with the one used in meson-exchange models. 

^We take the opportunity here to remark that the formula for F^^ reported in [9], though written for 
general warp factor a(z), is only correct in the case of AdS5 space in which a{z) = L/z. 
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0.68 fm 
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-1.461 


-1.459 


0.1% 



Table 1: Prediction of the nucleon observables with the microscopic parameters fixed by a 
fit on the mesonic observables. The deviation from the empirical data is computed using the 
expression \th — exp\/ mm{\th\, \exp\), where th and exp denote, respectively, the prediction 
of our model and the experimental result. 

Comparison with Experiments 

Let us now compare our results with real- world QCD, we therefore fix the number of colors 
Nc = 3 and choose our microscopic parameters to be 1/L ~ 343 MeV, M^L ^ 0.0165 and 
a ~ 0.94 (7 ^ 1.23). These values are obtained by minimizing the root mean square error 
(RMSE) in the mesonic sector. The detailed list of the observables we used can be found 
in [9] and the minimum RMSE for mesons is found to be 11%. 

The numerical results of our analysis and the deviation with respect to the experimental 
data are reported in table [Tj We find a fair agreement with the experiments, a 36% total 
RMSE which is compatible with the expected size of 1/Nc corrections. We discussed in the 
previous section that the isovector radii are divergent because of the chiral limit, it would 
be interesting to add the pion mass to the model and compute these observables. Table [T] 
also shows the proton-neutron magnetic moment ratio, which is in perfect agreement with 
the experimental value. Notice that for this observable, due to the different scalings of fis 
and /iy with Nc, our computation includes two orders of the l/N^ expansion: the leading 
order value which is —1 and the next-to-leading l/N^ correction which accounts for the extra 
—0.46. The axial charge is the one which shows the larger (almost 100%) deviation, and 
indeed removing this observable the RMSE decreases to 21%. We cannot exclude that, in 
a theory in which the naive expansion parameter is 1/3, enhanced 80% relative corrections 
to few observables might appear at the next-to-leading order. This failure in qa, therefore. 
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Figure 1: Scalar (left) and vector (right) electric form factors. We compare the results with 
the empirical dipole fit (dashed line) [7]. 




0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.0 0.1 0.2 0.3 0.4 0.5 0.6 

g- (GeV-) q- (GeV-) 

Figure 2: Normalized scalar (left) and vector (right) magnetic form factors. We compare the 
results with the empirical dipole fit (dashed line) [7]. 

does not invalidate the general picture. 

It is interesting to notice that a much better prediction for is obtained if one uses, 
instead of the standard procedure [6] considered in this paper, a different approach to the 
quantization of collective coordinates of the skyrmion, which has been proposed in Ref. [17]. 
The results of Ref. [17] can be directly applied to our case since, for what concerns the 
collective coordinate quantization, the 5D nature of our soliton is immaterial. We therefore 
find that the prediction for fis and for the radii are unaffected while both fiy and qa are 
rescaled by 5/3. We still obtain a good prediction for fxy = 2.98 (which is 27% away from 
the experimental value) and a much better prediction for qa = 1.17. Being the quantization 
of [17] equivalent to the standard one at large- A^^c, we have no reasons to prefer, a priori, one 
or the other. We have no reason either, however, to believe that the 1/Nc corrections one 
includes in this alternative approach really capture the leading 1/Nc corrections or at least 
part of them. If this was the case we should, of course, use the non-standard quantization 
and the discrepancy in the prediction of qa would disappear. 
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Figure 3: Left: deviation of the ratio of proton and neutron magnetic form factors from the 
large value (solid line), compared with the dipole fit of the experimental data (dashed 
line). Right: normalized axial form factor (solid line) compared with the empirical dipole fit 
(dashed lines) [7] and with the experimental data taken from [30,31]. 



If we stick, on the contrary, to the standard quantization procedure a small value of qa 
idA = 0.65 [6]) is also obtained in the original Skyrme model, but the situation improves if 
the effects of the p and u mesons are taken into account. The "complete" model described in 
Ref. [7] seems the one which should better mimic our 5D scenario, and qa = 0.99 in that case. 
The explicit chiral symmetry breaking, which is turned on in [7] , could explain the difference 
because the axial coupling is strongly sensitive to the large-r behavior of the solution (see 
the discussion following eq. (41 )) which is in turn heavily affected by the presence of the pion 
mass. Correction to qa from chiral symmetry breaking could therefore be enhanced. Notice 
that, however, this expectation fails in the original Skyrme model, where the addition of the 
pion mass does not affect qa significantly [27]. 

In figs, [g |2| and I we compare the normalized nucleon form factors at 7^ with the 
dipole fit of the experimental data. The shape of the scalar and axial form factors is of the 
dipole type, the discrepancy is mainly due to the error in the radii. The shape of vector form 
factors is of course not of the dipole type for small g^, but this is due to the divergence of the 
derivative at = 0. Including the pion mass will for sure improve the situation given that 
it will render finite the slope at zero momentum; it would be interesting to see if the dipole 
shape of these form factors is recovered in the presence of the pion mass. We also plot in the 
left panel of fig. |3]the deviation of ratio of the proton and neutron magnetic form factors from 
the large Nc value which is given, due to the the different large- A^^c scaling of the isoscalar 
and isovector components, by G^,j{q)/Gfj{q) = — 1. Not only we find that this quantity is 
quite well predicted, with an error < 15%, but we also see that its shape, in agreement with 
observations, is nearly constant away from = 0. Also in this case correction from the pion 
mass are expected to go in the right direction. 
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A The Equations of Motion 

In this technical appendix we report the EOM for the 2D fields which appear in our ansatz 
in eqs. ^ and (14) and we explain the notation used throughout the paper. 

The Residual Gauge Invariance 

Before discussing the detailed form of the EOM, it is useful to observe that our ansatz has 
not fixed the 5D gauge freedom completely, its form is indeed preserved by chiral SU{2)l,r 
gauge transformations of the form = U{t) ■ g ■ ?7^(t) and ql = U (t) ■ g'^ ■ ?7^(t) with 

g = exp[ia(r, z)xVa/(2r)] . (49) 



The operators A defined in eq. (10) have simple transformation rules under the residual 
symmetry. Indeed 

^^(i),afe^^/2^t = cosaA(i)''*Va/2 + sinaA(2).'^Va/2 , 
c/A(2)'«Va/2c/t = cosaA(2).«Va/2 - sinaA(i)'''V„/2 , 



so that the 2D fields (f)(x) and Xix) defined respectively in eq. ([9j) and (14) transform as charged 
complex scalars under this residual U (1). It is not hard to see that the fields Ap_ = {Ai, A2} 
transform as gauge field, so that under a residual transformation one has 

^ + dpa{r, z) , 
= 01 + e^"(^'^V , (50) 

while all the other fields are invariant. 
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There is also a second residual U{1) associated with chiral U{1)l,r 5D transformations 
of the form ^r = ^ and Tji = Tj'^ with 

{k ■ x) 



g = exp 



iP{r, z) 



a 



(51) 



Under this second residual U{1) only i?^ = {i?i,i?2} and p transform non trivially. We have 

and therefore -B^ is a gauge field and p a Goldstone. 

In order to make manifest the residual gauge invariance in the action and the EOM we 
introduced gauge covariant derivatives for the 0, x a^nd p fields 



r (2?p0)(.) = 9^0(.) + e(^^)A^0(,) 
( Dfip = dfip - Bp 



(53) 



The Equations of Motion 

The easiest way to derive the EOM for the 2D fields is to start from the Lagrangian and 
substitute the ansatz. Using the 5D action in eqs. ^ and ^ and rewriting the 5D fields 
in terms of the 2D ones (eqs. (11), (13) and ([li))), after a straightforward computation 



one finds the expressions for the mass M and moment of inertia A given in eqs. (16) and 



(17). Notice that in order to obtain the order terms of the action one has to perform a 
symmetric integration in d^x, which can simply be implemented by the replacement x^x^ -—>■ 



1/3(5*-'. We report here two contraction identities of the "doublet" operators A (eq. (10)) 
which can be useful for the computation of the 2D action 



^{x),ab^bid^^ _ ^(xy) ^{y),ai _ 



(54) 
(55) 



The EOM for the 2D fields can be simply obtained, at this point, by imposing the 
variation of the 2D action to vanish. We have also checked the consistency of our ansatz 
by showing that the same EOM are obtained by substituting directly into the 5D equations 



(12). The EOM for the fields which are already turned on in the static case are 



(a(z)Dn<h) + 



a{z) 



:i- 



{r'^a{z)A-^^) - a{z) {i(j)^D^(l) + h.c) + -iLe^^'d-^ 
df, {a{z)d~^s) - ^e'^^ [dp{-t(f)^D,(f) + h.c.) + A^p] = 



= 



(56) 
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while the equations for the "new" fields which are turned on for the rotating skyrmion are 



d''{r\{z)df,v) - 2a{z) [v{l + |0|') - x<P^ - <Px^] + 



l)Bf,^ + rQA 



fii/ 



D-^{r^a{z)Df,x) + a{z) [2t;0 - (1 + - 7 Le^' (D f^cj)) [td.irQ) + D,p] = 



7L 



Df,(f){D^xy + h.c.) + -Ai,,{2v - x<f)^ - 



dfi{a{z)Dfip) 



^Dfipdi, (- 
OL^ \r 

{D-,<P{D,x)^ + h.c) + Ja^p(0x^ - X0^) + \d-,{rQ)d, 



{r^a{z)Bp-^) + 2aiz)Dpp 

+7Le~^'{ [ix - v(P){D,^y + h.c] + (1 - \(P\^)d,v 



;Qdi 



(57) 



In order to solve numerically the EOM, they must be rewritten as a system of elliptic 
partial differential equations. This can be achieved by choosing a 2D Lorentz gauge condition 
for the residual U{1) gauge fields 



0. 



d'^Bn = . 



(5^ 



In this way the equations for Ap become J'^ = (r'^aAP-'^) = r'^adpd^'A'^ + dp{r'^a)A 
which is an elliptic equation and a similar result is obtained for B^. As discussed in [8], to 
impose the gauge condition, one can solve the "gauge-fixed" EOM counting the gauge field 
components as independent fields. In this way, if one imposes the gauge conditions at the 
boundaries, then the gauge is maintained also in the bulk. 



The Boundary Conditions 



The IR and UV boundary conditions on the 2D fields follow from eq. ^ and eq. ([s]) and 



from the gauge choice in eq. ( 58 ) . They are given explicitly by 



z — Zi^ 



and 



Z = Zuy 



01 = 

^202 = 

Ai = 
d^A^ = 

d2S = 

'01 = 

02 = -1 

Ai = 
d2A2 = 
s = 



Xi = 

d2X2 = 

d2V = 
d2Q = 






(59) 



(60) 
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The boundary conditions at r = oo have to ensure that the energy of the solution is 
finite, this means that the fields should approach a pure-gauge configuration. At the same 
time one has to require that the solution is non-trivial and its topological charge (eq. ([s])) is 
non zero. To obtain a soliton solution with B = 1 one can impose the conditions 



oo 






(61) 



The r = boundary of our domain requires an ad hoc treatment, given that the EOM 
become singular there. Of course this boundary is not a true boundary of our 5D space, 
but it represents some internal points. Thus we must require the 2D solution to give rise to 
regular 5D vector fields at r = and we must also require the gauge choice to be fulfilled. 
These conditions are 



(l + 02)/r- 
A2 = 
diAi = 
s = 





(62) 



B Numerical Techniques 



To obtain the numerical solution of the EOM we used the COMSOL 3.4 package [32], which 
permits to solve a generic system of differential elliptic equations by the finite elements 
method. A nice feature of this software is that it allows us to extend the domain up to 
boundaries where the EOM are singular [i.e. the r = line), because it does not use the 
bulk equations on the boundaries, but, instead, it imposes the boundary conditions. 

In order to improve the convergence of the program and the numerical accuracy, one is 
forced to perform a coordinate and a field redefinition. The former is needed to include the 
r = oo boundary in the domain in which the numerical solution is computed. The advantage 
of this procedure is the fact that in this way one can correctly enforce the right behaviour of 
the fields at infinity by imposing the r = oo boundary conditions. A convenient coordinate 
change is given by 

X = carctan ^-j , (63) 

where x is the new coordinate used in the program and c is an arbitrary constant. The 
domain in the x direction is now reduced to the interval [0,c7r/2]. The parameter c has 
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been introduced to improve the numerical convergence of the solution. A good choice for c 
is c ~ 10, which allows to have a reasonable domain for x and, at the same time, does not 
compress the solution towards a; = 0. 




A field redefinition is needed to impose the regularity conditions at r = (eq. (62)). For 
this purpose we use the rescaled fields 



(64) 



With these redefinitions, in the new coordinates, the r = boundary conditions read as 

'^rtr" f-Bi=o 

O: iXZl ^-j;: . ,65) 

= [ Q = 

In order to ensure the convergence of the program another modification is needed. As 
already discussed, to obtain a soliton solution with non-vanishing topological charge we have 



to impose non-trivial boundary conditions for the 2D fields at r = oo (eq. (61)). It turns 
out that, if such conditions are imposed, the program is not able to reach a regular solution. 
This is so because the r = oo boundary is singular and imposing non-trivial (though gauge- 
equivalent to the trivial ones) boundary conditions at a singular point spoils the regularity 
of the numerical solution; the same would happen if the topological twist was located at 
r = 0. To fix this problem we have to perform a gauge transformation which reduces the 
r = oo conditions to trivial ones and preserves the ones at r = at the cost of introducing 
a "twist" on the UV boundary. For this, we use a transformation of the residual U{1) chiral 



gauge symmetry associated to SU{2)L^ji (eq. (49)) with 



air,z) = il-z/L)fir), (66) 
where /(r) can be an arbitrary function which respects the conditions 

l^(°)=° and /nO) = m 

I /(oo) - vr [ r (oo) • ^^^^ 



For c ~ 10 it turns out that a good choice for /(r) is /(r) = 2arctanr. The gauge-fixing 

lified as 

drA, + d,A2 - (1 - z/L)f"{r) = , (68) 



condition for Ap, is now modified as 
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the UV boundary conditions are given by 



^ — ^^uv 



' X'?/'! = sin fir) 
(—1 + x'^'i) = — cos /(r) 
Ai = fir) 
dzA2 = 
s = 



and the r = oo constraints are now trivial 

r^i = o 



oo 



(-1 + X^2) 

d,Ai = 
^2 = 
s = 



Xi = -sin/(r) 
X2 = cos/(r) 
v = -l 
Q = 






(69) 



(70) 



whereas the r = and the IR boundary conditions are left unchanged. Notice that in the 



new gauge the EOM for Ap are modified in accord to eq. (68), however they are still in the 
form of elliptic equations. 
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